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Abstract 

We go over our finding tliat tlie dynamics at the noise-perturbed 
edge of cliaos in logistic maps is comparable to that observed in super- 
cooled liquids close to vitrification. That is, the three major features 
of glassy dynamics in structural glass formers, two-step relaxation, 
aging, and a relationship between relaxation time and configurational 
entropy, are displayed by orbits with vanishing Lyapunov exponent. 
The known properties in control-parameter space of the noise-induced 
bifurcation gap play a central role in determining the characteristics 
of dynamical relaxation at the chaos threshold. Time evolution is 
obtained from the Feigenbaum RG transformation, it is expressed an- 
alytically via g-exponentials, and described in terms of nonextensive 
statistics. 
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1 Introduction 

The problem of understanding the dynamics of glass formation is known to 
be difficult and challenging and therefore continues to be a frontier topic in 
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statistical physics [1]. It is generally thought that ordinary phase-space mix- 
ing is not wholly realized during glass forming dynamics, i.e. upon cooling, 
caged molecules rearrange so slowly that they cannot sample configurations 
in the available time allowed by the process [1]. The advance to the glass 
transition is signaled by very pronounced deceleration of relaxation processes 
since the characteristic times for their completion increase several orders of 
magnitude within a small interval of temperatures. Because of this extreme 
circumstances an important theoretical issue is to determine whether under 
conditions of ergodicity malfunction, and, as a final point, downright failure, 
the Boltzmann-Gibbs (BG) statistical mechanics is still capable of describing 
stationary states on the point of glass formation or those representing the 
glass itself. This question and also that about the possible applicability of 
generahzations of the BG statistics, such as the nonextensive statistics [2] 
[3], to glass formation is a subject of current concern. 

As a fresh approach we study the dynamics of a model in which only the 
essential ergodicity breakdown feature is present and where there is neither 
molecular structure nor molecular interactions. The general aim is to learn 
if the occurrence of the key properties of glassy dynamics are due primar- 
ily to the approach to a nonergodic state and therefore are observable in 
completely different types of systems, hence suggesting some degree of uni- 
versality. Our minimal model for glass dynamics is the prototypical logistic 
map at the edge of chaos but perturbed by small amplitude noise. It has 
only one degree of freedom but the consideration of external noise could be 
thought to represent the effect of many other systems coupled to it, like in 
the so-called coupled map lattices [4]. Recently [5], we have proved that 
the dynamics of this map shows robust similarities with that of supercooled 
liquids close to vitrification. This is so since two-step relaxation, aging, and 
a relationship between relaxation time and configurational entropy, are the 
three main (phenomenologically established) characteristics of glassy dynam- 
ics in structural glass formers [1]. The basic ingredient of ergodicity failure is 
obtained for orbits at the onset of chaos in the limit towards vanishing noise 
amplitude. Here expand the discussion of our results [5] . 

Our view is that glass formation is one in which the system is driven 
gradually into a nonergodic state by reducing its ability to pass through 
phase-space-filling configurational regions until it is only possible to go across 
a (multi) fractal subset of phase space. This situation is generated in the 
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logistic map with additive external noise, 

xt+i = fix{xt) = 1 - fJ'Xt + xtcr, -1 < < 1, < // < 2, (1) 

where Xt is Gaussian-distributed with average {xtXf) = ^t.t', and a measures 
the noise intensity [12] [13]. Notice that the formula (1) can be written as a 
discrete form for a Langevin equation. 

2 Dynamics of glass formation 

Wc recall the main dynamical properties displayed by supercooled liquids on 
approach to glass formation. One is the growth of a plateau and for that 
reason a two-step process of relaxation, as presented by the time evolution 
of correlations e.g. the intermediate scattering function Fk(At) [1]. This 
consists of a primary power-law decay in time difference At (so-called f3 
relaxation) that leads into the plateau, the duration of which diverges 
also as a power law of the difference T — Tg as the temperature T decreases 
to a glass temperature Tg. An 'ideal' glass transition has been conjectured 
to occur at a Tc which is not experimentally accessible but that approaches 
Tg at infinitely slow cooling in fragile glass formers [14]. After there is a 
secondary power law decay (so-called a relaxation) away from the plateau 
[1]. A second important (nonequilibrium) dynamic property of glasses is the 
loss of time translation invariance observed for T below Tg, a characteristic 
known as aging [6] , that is due to the fact that properties of glasses depend 
on the procedure by which they are obtained. Remarkably, the time fall off 
of relaxation functions and correlations display a scaling dependence on the 
ratio t/tyu where is a waiting time. A third notable property is that 
the experimentally observed relaxation behavior of supercooled liquids is 
effectively described, via reasonable heat capacity assumptions [1], by the 
so-called Adam-Gibbs equation, = Aexp{B /TSc), where the relaxation 
time can be identified with the viscosity or the inverse of the difussivity, 
and the configurational entropy Sc is related to the number of minima of 
the fluid's potential energy surface (and A and B are constants) [1] . A first- 
principles derivation of this equation has not been developed at present. As 
the counterpart to the Adam-Gibbs formula, we show below that our one- 
dimensional map model for glassy dynamics exhibits a relationship between 
the plateau duration tx, and the entropy Sc associated to the iterate positions 
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(configurations) within the largest number of phase-space bands allowed by 
the bifurcation gap - the noise-induced cutoff in the period-doubling cascade 
[12]. 

Our results suggest that the properties known to be basic of glassy dy- 
namics [1] are likely to manifest too in completely different physical problems, 
such as in nonlinear dynamics, e.g. the mentioned coupled map lattices [7], in 
critical dynamics [8] , and other fields. This hints that new predictions might 
be encountered in the studies, experimental or otherwise, of slow dynamics 
displayed by systems other than liquids close to the glass transition. In re- 
lation to this, aspects of glassy dynamics have been observed in metastable 
quasistationary states in microcanonical Hamiltonian systems of classi- 
cal rotors with homogeneous long-ranged interactions. For special types of 
initial conditions it has been found that both two-step relaxation [9], where 
the length of the metastable plateau diverges with infinite size N ^ oo, and 
aging [10], [11] are present in these systems. In our simpler one-dimensional 
dissipative map the amplitude a plays a role parallel to T — Tg or T — Tc 
in the supercooled liquid or 1/N in the system of rotors. Notice that the 
equivalence between a and T — Ty is not literal as a cannot take negative 
values. 

3 Dynamics at the edge of chaos 

In the absence of noise cr = the Feigenbaum attractor at /x = yUc(O) = 
1.40115... is the accumulation point of both the period doubling and the 
chaotic band splitting sequences of transitions [12]. Except for a set of zero 
measure, all the trajectories with /ic(0) and initial condition —1 < Xin < 1 fall 
into the attractor with fractal dimension df — 0.5338.... These trajectories 
represent nonergodic states, since as t ^ oo only a Cantor set of positions is 
accessible out of the total phase space — 1 < a; < 1. For cr > the noise fiuc- 
tuations wipe the sharp features of the periodic attractors as these widen into 
bands similar to those in the chaotic attractors, nevertheless there remains 
a well-defined transition to chaos at yUc(o") where the Lyapunov exponent Ai 
changes sign. The period doubling of bands ends at a finite value 2^*^*^^ as the 
edge of chaos transition is approached and then decreases at the other side of 
the transition. This effect displays scaling features and is referred to as the 
bifurcation gap [12] [13]. When a > the trajectories visit sequentially a set 
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of 2" disjoint bands or segments leading to a cycle, but the behavior inside 
each band is completely chaotic. These trajectories represent ergodic states 
as the accessible positions have a fractal dimension equal to the dimension 
of phase space. Thus the removal of the noise cr — > leads to an ergodic to 
nonergodic transition in the map and we compare its properties with those 
known for the process of vitrification of a liquid as T — T^. 

The manner in which the attractor at the onset of chaos jJ^dO) is visited in 
time was analyzed recently with the use of the initial condition Xj„ — [15]. 
It was found that the absolute values for the positions Xr of this trajectory at 
time-shifted r = t + 1 has an structure consisting of subsequences with a com- 
mon power- law decay of the form r^'^/'^^i with g=l — ln2/lnQ;~ 0.24449, 
where a = 2.50290... is the Feigenbaum universal constant that measures the 
period-doubling amplification of iterate positions [15]. That is, the Feigen- 
baum attractor can be decomposed into position subsequences generated 
by the time subsequences r = {2k + 1)2"', each obtained by proceeding 
through n = 0,1,2,... for a fixed value of A; = 0,1,2,.... The A; = sub- 
sequence can be written as Xt = exp2_g(— Agt) with Ag = In a/ In 2, and 
where expg{x) = [1 — {q — l)x]^^^~^ is the g-cxponential function. These 
properties follow from the use of Xin = in the scaling relation [15] 



g'^'ixir^ =^-1/1-5 .r.i/W. 



^(r^/^-%.„) . (2) 



In the presence of noise (u small) one obtains instead [5] 

g{r'/'-^x) + xar'/'-'GAir'/'-'x) , (3) 



where G\{x) is the first order perturbation eigenf unction, and where r = 
1 - ln2/lnK ~ 0.6332. Use of Xi^ = yields Xr = r-^i-? |i + ^^o-rVi-'"! or 
Xt — exp2-q{—\t) [1 -|- X(jexp^(Art)] where t — T — 1 and A^ = In In 2. 

At each noise level a there is a 'crossover' or 'relaxation' time t^, = r^, — 1 
when the fiuctuations start suppressing the fine structure of the orbits with 
Xin = 0. This time is given by Tx = cr^'^, the time when the fiuctuation 
term in the perturbation expression for Xr becomes unbounded by a, i.e. 
Xr^ — T~^/^~'^ |1 There are two regimes for time evolution at ndcr). 

When T < Tx the fiuctuations are smaller than the distances between neigh- 
boring subsequence positions of the a = orbit at /Uc(0), and the iterate 
positions with a > fall within small non overlapping bands each around 
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the cr = position for that r. Time evolution follows a subsequence pat- 
tern analogous to that in the noiseless case. When t Tx the width of the 
noise-generated band reached at time — 2^ matches the distance between 
adjacent positions where N ~ — Incr/ln/t, and this implies a cutoff in the 
progress along the position subsequences. At longer times r > r^. the orbits 
no longer follow the detailed period-doubling structure of the attractor. The 
iterates now trail through increasingly chaotic trajectories as bands merge 
with time. This is the dynamical image - observed along the time evolution 
for the orbits of a single state /Uc(o') - of the static bifurcation gap originally 
described in terms of the variation of the control parameter [13], [16], [17]. 
The plateau structure of relaxation and the crossover time can be clearly 
observed in Fig. lb in Ref. [18] where < > — < Xt is shown for several 
values of a. 



4 Parallels with glassy dynamics 

At noise level a the orbits visit points within the set of 2^ bands and, as 
explained in Ref. [5], this takes place in time in the same way that period 
doubling and band merging proceeds in the presence of a bifurcation gap 
when the control parameter is run through the interval < /x < 2. Namely, 
the trajectories starting at Xin — duplicate the number of visited bands at 
times r = 2", n = 1, A^, the bifurcation gap is reached at Tx — 2^, after 
which the orbits fall within bands that merge by pairs at times r = 2^+", 
n = 1,...,N. The sensitivity to initial conditions grows as = exp^^Xgt) 
[q = 1 — In 2/ In q; < 1) for t < tx, but for t > tx the fluctuations dominate 
and grows exponentially as the trajectory has become chaotic {q — 1) 
[5]. This behavior was interpreted [5] to be the dynamical system analog 
of the a relaxation in supercooled fluids. The plateau duration t^. — cxd as 
cr — > 0. Additionally, trajectories with initial conditions Xin not belonging to 
the attractor exhibit an initial relaxation process towards the plateau as the 
orbit approaches the attractor. This is the map analog of the P relaxation 
in supercooled liquids. 

The entropy Sc{fic{<^)) associated to the distribution of iterate positions 
(configurations) within the set of 2^ bands was determined in Ref. [5]. 
This entropy has the form 5'c(/Xc(c)) — 2^as, since each of the 2-^ bands 
contributes with an entropy as, where s — — X^iP(x) lnp(x)dx and where 
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p{x) is the distribution for the noise random variable. Given that 2^ = l+t^ 
and a = (1 + tx)^^^^^'^, one has ScifJ^ctx)/ s = (1 + ta;)^'^''"'^^'''or, conversely, 

tx = {s/S^Y'-^y^. (4) 

Since tx cr''"^, r — 1 ~ —0.3668 and (1 — r)/r 0.5792 then tx ^ oo and 
— > as (7 — > 0, i.e. the relaxation time diverges as the 'landscape' entropy 
vanishes. We interpret this relationship between tx and the entropy Sc to be 
the dynamical system analog of the Adam-Gibbs formula for a supercooled 
liquid. Notice that Eq.(4) is a power law in while for structural glasses 
it is an exponential in [1]. This difference is significant as it indicates 
how the superposition of molecular structure and dynamics upon the bare 
ergodicity breakdown phenomenon built in the map modifies the vitrification 
properties. 

The aging scaling property of the trajectories Xt at /Uc(a") was examined 
in Ref. [5]. The case o" = is readily understood because this property is 
actually built into the position subsequences Xt = g^'^\0) , r = {2k + 1)2", 
k,n = 0,1, ... referred to above. These subsequences can be employed for the 
description of trajectories that are at first held at a given attractor position 
for a waiting period of time and then released to the normal iterative 
procedure. For illustrative purposes we select the holding positions to be 
any of those for a waiting time = 2k + 1, k = 0,1, ... and notice that for 
the — orbit these positions are visited at odd iteration times. The lower- 
bound positions for these trajectories are given by those of the subsequences 
at times {2k + 1)2". See Fig. 1 in Ref. [5]. Writing r as r — tu, + t we have 
that t/ty, = 2" - 1 and Xt+t^ = ^(*-)(0)c/(*/*™)(0) or 

Xt+t^^9^'-\0)e^g{-X,t/t^). (5) 

This fully developed aging property is gradually modified when noise is 
turned on. The presence of a bifurcation gap limits its range of validity to 
total times tw + t < tx{(j) and so progressively disappears as a is increased. 
Recently [19] aging has been observed in the subdiffusion associated to a 
nonlinear map that, as it is the case here, has Ai = [20]. 

5 Discussion 

Thus, the dynamics of noise-perturbed logistic maps at the chaos threshold 
exhibits the most prominent features of glassy dynamics in supercooled liq- 
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uids. The existence of this analogy cannot be considered accidental since 
the limit of vanishing noise amplitude a — > (the counterpart of the limit 
T — — > in the supercooled hquid) entails loss of ergodicity. The incidence 
of these properties in such simple dynamical systems, with only a few degrees 
of freedom and no reference to molecular interactions, suggests a universal 
mechanism underlying the dynamics of glass formation. As definitely proved 
[15], the dynamics of deterministic unimodal maps at the edge of chaos is 
a genuine example of the pertinence of nonextensive statistics in describing 
states with vanishing ordinary Lyapunov exponent Ai. Here we have shown 
that this nonergodic state corresponds to the limiting state, a ^ 0, ^ oo, 
for a family of small a noisy states with glassy properties, that are noticeably 
described for t < tx via the g-exponentials of the nonextensive formalism [15]. 

To ascertain the degree of parallelism between the map at /Xc(cr) and a 
thermal system we keep in mind the following criteria: The exponential 
(or Ai > 0) of a chaotic state is the counterpart of the BG equilibrium state. 
On the other hand, a power-law (or Ai = 0) of an incipient chaotic state 
represents an "anomalous stationary state", q 7^ I, and in the case of /ic(0) 
a nonergodic state. Only when cr = we have a true anomalous stationary 
state (of infinite duration) that exhibits a precise aging property. In the map 
at /ic(o") we find two regimes separated by with presence of aging when 
t < tx- In glass formers the two-step relaxation is seen in equilibrium two- 
time correlations while aging is a manifestation of the system having fallen 
out of equilibrium. Because of these apparent differences a more detailed 
comparison between the map and the thermal system properties is required, 
and it may be necessary to consider neighbouring states around /U^(o"). These 
states would be expected to display also slow dynamics. When /i > /Uc((j) 
one has Ai > and the analog of BG equilibrium, whereas for /i < ndcr) one 
has Ai < and a situation "out of equilibrium". 

The anomalous stationary state has properties akin to both equilibrium 
and out of equilibrium. Equilibrium-like since it is stationary, an orbit that 
starts in the attractor remains there forever. Noncquilibrium-like because 
the finite-time Lyapunov exponents fluctuate and take negative values (im- 
plying confinement in phase space). When o" > it would be expected that 
the frequency of these negative values decreases sharply after t^- Further 
studies arc needed here. Also, there remains the question about whether 
some features of glassy dynamics, such as loss of time translation invariance, 
are due to nonergodicity or to merely being out of equilibrium (i.e. in a 
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non-stationary state). 

It has been suggested on several occasions [3] tliat the setting in which 
nonextensive statistics appears to emerge is hnked to the incidence of nonuni- 
form convergence, such as that involving the thermodynamic N ^ oo and 
very large time t ^ oo limits. For example, in the rotor problem mentioned 
above - for specific choices of initial conditions - if iV — oo is taken before 
t — > oo the anomalous metastable states with noncanonical properties appear 
to be the only observable stationary states, whereas if t ^ oo is taken before 
N ^ oo the usual BG equilibrium states are obtained. Here it is clear that 
a similar situation takes place, that is, if c ^ is taken before t — oo a 
nonergodic orbit confined to the Feigenbaum attractor is obtained, whereas 
if t — > cxD is taken before cr — > a typical q = I chaotic (ergodic) orbit is 
observed. 

Finally, it is worth mentioning that while the properties displayed by the 
map capture in a qualitative, heuristic, way the phenomenological issues of 
vitrification, they are obtained in a quantitative and rigorous manner as the 
map is concerned. Our map setup is a rarely available 'laboratory' where 
every aspect of the dynamics can be studied analytically. 
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